Suppression of friction by mechanical vibrations 
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Mechanical vibrations are known to alfect frictional sliding and the associated stick-slip patterns causing 
sometimes a drastic reduction of the friction force. This issue is relevant for applications in nanotribology and to 
understand earthquake triggering by small dynamic perturbations . We study the dynamics of repulsive particles 
confined between a horizontally driven top plate and a vertically oscillating bottom plate. Our numerical results 
show a suppression of the high dissipative stick-slip regime in a well defined range of frequencies that depends 
on the vibrating amplitude, the normal applied load, the system inertia and the damping constant. We propose 
a theoretical explanation of the numerical results and derive a phase diagram indicating the region of parameter 
space where friction is suppressed. Our results allow to define better strategies for the mechanical control of 
friction. 

PACS numbers: 81.40.Pq, 46.55.+d, 68.35.Af, 68.08.-p 



Natural or artificially induced manipulations by small me- 
chanical vibrations, when applied at suitable frequency and 
amplitude ranges, may help in driving a contacting sliding 
interface out of its potential energy minima, thus increas- 
ing considerably surface mobility and diffusion, and reducing 
friction. This has been shown experimentally for sliding of 
nanoscale contacts through, e.g., the atomic force microscope 
Qfllstl. ™d in computer simulations via extended molecu- 
lar dynamics [4J and simple modeling approaches [5, 6^^. 
On a larger scale, it has been observed that in sheared gran- 
ular media experiments the stick slip behavior is significantly 
perturbed by tiny transverse vibrations |&^9J. Since geologi- 
cal faults are often filled with a granular gouge, these results 
might be relevant to understand earthquake triggering by low 
amplitude seismic waves [lOJ. Despite these promising nu- 
merical and experimental contributions, a quantitative theory 
accounting for the friction dependence on vibrations is still 
lacking. 

In this letter, we study the frictional properties of a two 
dimensional set of repulsive particles confined between two 
rigid plates (see FigdJ 12]. The top plate is attached to 
a spring that is pulled at constant velocity, while the bottom 
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Oscillating bottom substrate Zbot = A sin(Lc;oi) 

FIG. 1 : Sketch of the system, with the rigid top and bottom plates in- 
dicated in black and the confined particles in red. The top is dragged 
through a spring of elastic constant K moving at constant velocity 
Vext, while the bottom vibrates vertically with frequency luo and am- 
plitude A 



plate is vibrated vertically. Without vibration, the top plate 
would slide exhibiting a characteristic stick-slip behavior Vi- 
brations induce a drastic reduction of the friction coefficient 
and a suppression of the stick-slip behavior, but only in a well 
defined frequency range. We propose a theoretical argument 
to explain this behavior and construct a phase diagram indi- 
cating the parameter region for which friction is suppressed. 
The theoretical results are in excellent quantitative agreement 
with numerical simulations. Finally, we investigate the fric- 
tional response of the system to a short vibration pulse and 
find large slip events when the vibration frequency lies in the 
appropriate regime. These results could be relevant to under- 
stand triggering of frictional slip by small perturbations. 

Our two dimensional system consists of two identical top 
and bottom substrates, composed of fit — rn, particles with 
coordinates r* 
aration 0^=1. 

between the top and bottom plates. The mass m of all the par- 
ticles is the same, so that the total mass of the confined layers 

is Mr, — m 



and rj' respectively and constant lattice sep- 
We confine Up particles with coordinates 
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while the top and bottom plate have a mass 
- m ■ Uf. Due to the vertical confinement, 
periodic boundary conditions are applied only along the x di- 
rection. The particles interact via a pairwise repulsive poten- 
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for r < ro and U{r) = 



otherwise. The parameters of the potential are the same for 
all the particles, wether in the top and bottom plate or in be- 
tween them. We adjust the number of particles rip, imposing 
that Up — Ni ■ rit — Ni ■ rib with Ni an integer number. In this 
way the systems forms Ni perfectly ordered layers of confined 
particles. We also consider the case of a few missing particles, 
finding exactly the same results as long as the system remains 
ordered. The top plate of mass is subject to a normal force Fm 
and is pulled through a spring, of elastic constant K, which 
moves along the horizontal x direction with constant velocity 

Vext 
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{Xhot, Ztot) the center of mass coordinates of the top and bot- 

torn plate respectively, where Utop = J2 ^1' ^> 

i=i 

rib 

— ^ rj, the particles satisfy the equations of motion 
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where Np + rif + rib and 77 is the damping coefficient that 
accounts for a viscous dissipation. The temperature is con- 
trolled by a Langevin thermostat according to the relation 
(f (t)'''^"f (t')™") = ^mr]kBTS{t - t'). In the present sim- 
ulations, we consider very low temperatures fc^T = IO^^C/q- 

To study the influence of mechanical vibrations on the sys- 
tem, the bottom plate is vibrated vertically Zbot = ^0 + 
Asin(woi) where Zq is a reference coordinate, A is the am- 
plitude and ojQ the frequency, while its horizontal component 
Xbot is held fixed. In all the simulations, we compute the in- 
stantaneous friction force Fl by measuring the spring elon- 
gation of the driving apparatus, {Xtop{t) — Vextt), so that 
Fl — K {Xtop{t) — Vextt)- The friction coefficient is defined 
as /Lt = Fi/ Fn and its average (ft) , obtained integrating its 
value over a sufficiently long time interval in the steady state. 

In Fig. 12^, we show the behavior of the friction coefficient 
/J, as a function of time in the case of three confined parti- 
cle layers for three distinct values of the frequency ujq. For 
ujQ — 1 and llIq = 3, we observe a characteristic stick slip 
behavior, with loading phases where the system is stuck, fol- 
lowed by rapid slip events in which the force accumulated by 
the spring is relaxed. A very similar pattern also takes place 
in absence of vibrations at low external driving rates, with the 
system alternately sticking and slipping forward (not shown). 
At ujQ = 2, we observe a strong reduction of the friction co- 
efficient and a drastic suppression of the sawtooth stick-slip 
behavior To make this observation more quantitative, we re- 
port in Fig. ^ and 12}; the systematic variations of the aver- 
age friction coefficient (fi) with the vibration frequency ujq. 
The left panel shows results obtained for various oscillation 
amplitudes A, ranging from 3% to 9% of the film thickness, 
while the right panel displays results for different values of the 
damping coefficient r/. We see that the suppression of friction 




FIG. 2: a) Friction coefficient fi vs time obtained for three values 
of the bottom oscillating frequency uiq (three lubricant layers), b) 
Time-averaged value of /i as a function of luq for four values of the 
oscillation amplitude A (ranging from 3% to 9% of the Ni = 3 
lubricant thickness), c) Time-averaged value of jj, vs ujo for Ni — 4 
and four values of the damping coefficient rj. 



appears in a well defined range of frequencies [wi, 0^2], which 
depends on A and rj. 

To understand the origin of this phenomenon, we compute 
the power spectrum of the vertical position of the top sub- 
strate Ztop for the three values of ojq corresponding to Fig|2^. 
At low frequencies (inset in Figl3t), the top plate and the con- 
fined particles vibrate in phase with the oscillations of the bot- 
tom plate. Hence, the spectrum displays a peak at frequency 
luq and its integer multiples and vibrations have no visible ef- 
fect on sliding friction. The situation is drastically different 
at cjo = 2 (Fig[3j>)- In this case, the top plate and the con- 
fined particles can not follow the bottom plate. The vertical 
position of the top plate increases, presenting high amplitude 
oscillations which diminish the contact time between the con- 
fined particles and the bottom plate, reducing considerably the 
friction force. The corresponding Fourier spectrum shows ad- 
ditional peaks at integer multiples of ujq /2. We find that this is 
a general feature of the spectrum in the interval of friction sup- 
pression [oji , W2]. Additional peaks may appear at frequencies 
uj — luq/uq, where hq is an integer, implying that the top plate 
oscillates with a period that is an integer multiple of the driv- 
ing period on the bottom plate. Further increases of the oscil- 
lation frequency induce a reduction of the ampUtude of the top 
plate oscillations around the equilibrium position (see FigO};), 
so that the low frequency noise suppresses the resonant peaks. 
As a result, the friction force increases again. 

The mechanism illustrated by the spectral analysis suggests 
an argument to derive analytically the values of uji and uj2, 
defining the interval of friction suppression. At low frequen- 
cies, the particles follow coherently the oscillation of the bot- 
tom plate owing to the downward action of the normal load 
Fn. The reduction of friction occurs when the particles are 
able to detach from the bottom plate due to the action of the 
inertial force induced by the external vibration Fin = MZ^ot, 
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FIG. 3: Power spectrum of the vertical component Zto-p of the plate 
for different values of the oscillating frequency ojo and A'^; = 3, 
A = 0.3, T] = \. At low frequency (ljo=1) the system follows 
the oscillating bottom substrate. At higher frequencies loq — 2 the 
system oscillates also at one-half frequency of uiq. At loq — 3 a low 
frequency noise appears. In the insets we report the time evolution 
of Ztop (dotted line) and Ztot (thick line). 



where AI = M.„ 
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The frequency uji corresponds to 



the condition for which the inertial force overcomes the com- 
bined action of the normal load and the damping force 
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Estimating the inertial force as Fin — MAoJq and the damp- 
ing force as F^amp — MpTjAajQ, we obtain an implicit expres- 
sion for the starting frequency of friction suppression uji 



MAluI = Fn + MpTjAuji. 
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It is convenient to work with dimensionless quantitities, defin- 
ing 
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Using these rescaled variables, Eq.|5]yields: 
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To estimate the recovery frequency u!2, we determine the 
conditions for the presence of low frequency vibrations. Due 
to the external oscillations, the confined layers detach from 
the bottom substrate during a characteristic time At that can 
be estimated as 



At ~ ZbotM/FN ^ AlooM/Fn 



(8) 



To observe low frequency vibrations, which will cause the re- 
covery of the friction force, the period of the external oscilla- 
tion should be smaller than the rise time associated with the 
internal vibrations of the particles (i.e. ^ < At). This condi- 
tion corresponds to the maximum of the momentum transfer 
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FIG. 4: a) Comparison between the numerical results (symbols) and 
the theory (solid line). u)\ as a function of the dimensionless vari- 
ables / and m. b) lu2 as a function of /. The behaviors correspond- 
ing to different symbols are obtained keeping three parameters fixed 
and varying the fourth, c) Phase diagram indicating the region of 
friction suppression. Top curve: UJ2 vs /. Grey-striped region: lo-l vs 
/ for I < m < 1. Friction is strongly suppressed for Cu-i < uj < Lb2 
(grey region). 



from the vibrating plate to the confined particles. Using again 
the dimensionless variables defined in Eq. (|6]l, we estimate 

(9) 

The theoretical predictions for loi and L02 are in excellent 
agreement with the numerical simulations, as shown in Fig. 

and |4j5. The numerical values are obtained varying the 
number of layers Ni, the vibration amplitude A, the damping 
coefficient 77 and the normal load Fm- Notice that the theory 
has no adjustable parameters. 

From the analytical relations (|7]i and (|9]l, we can draw a 
phase diagram indicating, in the space of dimensionless vari- 
ables (|6l), the region where friction is suppressed (Fig|4};). The 
region of friction reduction is enclosed between uJi and Cj2 
and shrinks as we reduce /, until it finally disappears. Notice 
that UJl depends on the value of reduced mass which lies be- 
tween rh = 1/2, for a single confined layer, and rh ~ \ for 
an infinitely wide system. The possible values for the oji are 
indicated by the grey-striped region in Fig|4]:. 

The phase diagram displayed in Fig|4j; indicates the fric- 
tional behavior of the system under a steady vertical oscil- 
lation. There are cases, however, where the external vibra- 
tion acts only for a small amount of time, such in the case 
of earthquakes or avalanches triggered by seismic waves. To 
address this issue, we analyze changes in the stick-slip pat- 
tern for small vertical vibration of finite duration Ty and fre- 
quency Lo, with Ty 3> l/o-io. To avoid discontinuities, we 
switch the perturbation on the bottom plate smoothly: Z^ot = 
Zq + f{t,Ty)Asa\{uJot), where f{t,Ty) = (tanh(t/r) - 
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media and nanoscale interfaces. Further work in this direc- 
tion could be useful to optimize friction control in technolog- 
ical nanodevices and to design better strategies to forecast the 
triggering of instabilities in materials and geosystems. 
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FIG. 5: The effect of a short time vibration pulse on the stick-slip 
behavior for different frequencies uuo, Ni = 3, A — 0.2. The shaded 
light grey regions indicate the duration of vibration pulse. The largest 
effect is visibile when the frequency lies in the range of uji < ujq < 

UJ2- 

tanh((< - Ty)/T))/2, with t < T„. If we chose Ty to be 
of the same order of magnitude as the stick time, we observe 
that the perturbation typically leads to small changes in the 
slip patterns. When oji < ojq < uj2, however, the systems 
exhibit a large slip event (see Fig.|5]). After the perturbation is 
removed the system recovers the original stick-slip behavior, 
without any long-range memory effects. Our results suggest 
that catastrophic events are more likely to be triggered when 
the perturbation lies in a definite frequency interval. 

In conclusions, we have clarified the role of vibrations in 
the frictional sliding of a confined system. The general mech- 
anism for friction suppression that we have uncovered is based 
on the reduction of the effective interface contacts produced 
by vibrations. Since the results depend only on the relation 
between inertial and dissipative forces, we expect them to be 
valid for a wide class of sliding systems, including granular 
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